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Abstract 

The probabilities for gaps in the eigenvalue spectrum of the finite dimension N x N 
random matrix Hermite and Jacobi unitary ensembles on some single and disconnected 
double intervals arc found. These are cases where a reflection symmetry exists and the 
probability factors into two other related probabilities, defined on single intervals. Our 
investigation uses the system of partial differential equations arising from the Fredholm 
determinant expression for the gap probability and the differential-recurrence equations 
satisfied by Hermite and Jacobi orthogonal polynomials. In our study we find second 
and third order nonlinear ordinary differential equations defining the probabilities in the 
general N case. For iV = 1 and iV = 2 the probabilities and thus the solution of the 
equations are given explicitly. An asymptotic expansion for large gap size is obtained 
from the equation in the Hermite case, and also studied is the scaling at the edge of the 
Hermite spectrum as ^ oo, and the Jacobi to Hermite limit; these last two studies 
make correspondence to other cases reported here or known previously. Moreover, the 
differential equation arising in the Hermite ensemble is solved in terms of an explicit ra- 
tional function of a Painleve-V transcendent and its derivative, and an analogous solution 
is provided in the two Jacobi cases but this time involving a Painleve-VI transcendent. 



AMS classification - primary: 15A52, secondary: 34A34, 34A05, 33C45 
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1 Introduction 



It is a celebrated discovery of Jimbo, Miwa, Mori and Sato |14] that the probabihty of an 
eigenvalue free region in the bulk of the infinite GUE (Gaussian unitary ensemble of random 
Hermitian matrices) can be expressed exactly in terms of a Painleve-V transcendent. This has 
initiated a number of studies ||l^, ||, 11] which, in addition to clarifying the general setting of 
the exact result, give formalisms that allow analogous results to be obtained in other cases. 
For example, the probability of an eigenvalue free region at the edge of the infinite GUE 
(appropriately scaled) has been expressed in terms of the Painleve-II transcendent |jl^ , while 
a Painleve-III transcendent has been shown to determine the distribution of the smallest 
eigenvalue in the infinite LUE (Laguerre unitary ensemble of non-negative matrices of the 
form A^A with A complex |[l^]). Moreover, for the finite classical ensembles with unitary 
symmetry, the probability of a single eigenvalue free region which includes an endpoint of the 
support of the weight has been expressed in terms of solutions of certain non-linear equations 



[18|. 



It is the objective of this work to provide three new evaluations of gap probabilities for 
particular finite classical random matrix ensembles with unitary symmetry. The ensembles 
considered are the finite GUE and the symmetric Jacobi unitary ensemble (JUE). We recall 
the eigenvalue p.d.f. for an ensemble with unitary symmetry is of the form 

N 

iiw2{xi) n (1.1) 

1=1 l<j<k<N 

where the weight function W2{X) determines the specific unitary ensemble: 

, e~^', GUE 
W2{X) = { ^ (1.2) 

(1-A)°(1 + A)^, JUE 

(the symmetric Jacobi ensemble refers to the case a = /? of the JUE). These ensembles can be 
realised in terms of matrices with independent Gaussian elements. For the GUE, each matrix 
X say must be Hermitian and have its diagonal elements xjj (which must be real) and upper 
triangular elements Xjk = ujk + ivjk chosen with p.d.f. N(0, l/\/2) and N(0, 1/2) + iN(0, 1/2). 
For the JUE one first constructs Q auxihary rectangular Mi x N and M2 x A'" (Mi, M2 > A'") 
matrices a and b respectively, with complex elements independently distributed according to 
N(0, l/\/2)-hiN(0, 1/V2). Then with A = a'^a, B = b^b, it can be shown that the distribution 
of the eigenvalues of A{A + B)~^ is an example of (|l.l| ) with 

u,2(A) = A^^i-^(1- A)^^2"^, 0<A<1. (1.3) 

The change of variables A 1-^ (1 — A)/2 shows that this is an example of the JUE. 

With £'2(0; /; tL'2(A); A^) denoting the probability that there are no eigenvalues in the 



interval / of an ensemble with eigenvalue p.d.f. (1.1), the probabilities to be calculated are 

S2(0;(-oo,-s) U (s,oo);e-^';iV) (1.4) 

and 



E2iO;{-l,-s)U{s,l);{l-X^r;N), E2{0;{-s,sy,{l - X^r^N). (1.5) 



The quantity ( |1.4| ) gives the probability that there are no eigenvalues in the GUE with 
modulus greater than s, as does the first quantity in (|l.5| ) for the JUE. The final quantity in 



2 



(|1.5D gives the probability that there are no eigenvalues in the JUE with modulus less than 
s (the analogous probability for the GUE has previously been evaluated 

These probabilities are special because, like the situation in which / is a single interval 
which includes an endpoint of the support of the weight function noted above, we will show 
that they can be evaluated exactly in terms of the solution of certain non-linear equations. 
Another motivating factor for our study is a recently discovered ||9| general identity satisfied 
by the probability i?2(0; /; W2{X); N) applicable whenever W2{\) is even and I is symmetrical 
about the origin. The identity states 

E2{0; I; W2{Xy, N) =E2{0; 1+ ;y-^/^W2{y^'^); [{N + l)/2]) 

x£;2(0;/+;y+i/2«;2(y'/');m) (1.6) 

where denotes the portion of / on R"*" in the variable y = and W2{y^^'^) '.= for y < 0. 
The probabilities (|1.4|) and (|1.5D are all of the form required by the LHS of this identity, and 
are thereby related to 

E2{0;{s\ocy,y^'/^e-y;m) (1.7) 

and 

E2{0; (s^ 1); y^i/2(i _ y^a. ^^(g. ^2). ^^1/2(1 _ y-^a. ^) (18) 

(m = [{N + l)/2], [A^/2]) respectively. The weight functions in ( |1.7| ) and ( |1.8| ) are again 
classical, but now the eigenvalues are excluded from a single interval which includes an 
endpoint of the support of the weight function. The corresponding probabilities are known 



in terms of Painleve-V and VI transcendents from the works of and |11]. Thus as a 
by-product of providing independent evaluations of the LHS of (|1.6D, we will be providing 
inter-relationships between solutions of nonlinear equations. This theme is further developed 
in m. 



In Section 2 the formalism of Tracy and Widom |18] giving coupled differential equations 
for the gap probability and some auxiliary quantities is revised. In Section 3 these coupled 
equations are reduced to a single ordinary differential equation specifying the probability 

(1.4) , and this equation is used to compute the large gap size behaviour of the probability, as 
well as various scaling limits. A similar study is undertaken in Section 4 for the probabilities 

(1.5) . In Section 5 the solution of the second order equations of Section 3 (Hermite case) and 
third order equations of Section 4 (Jacobi case) are given in terms of certain Painleve-V and 
Painleve-VI transcendents respectively. 
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2 The General Formalism 



We are going to consider a set of the eigenvalue spectrum arising from the matrix ensembles 
consisting of an arbitrary number of disconnected intervals. In general the eigenvalues may 
be excluded from M disjoint intervals which together form /. Thus with the endpoints of 
these intervals denoted {aj}^^^, 

M 

1= IJ (a2m-l,a2m) • (2.1) 
m>l 

The probability of no eigenvalues being found in this interval is given by the general expression 
(see e.g. PI) 



n=l ■ 



dxi... J dXnPn{xi, ... ,Xn) , (2.2) 

where pn is the n-point distribution function of the eigenvalue p.d.f. For matrix ensembles 
with unitary symmetry the eigenvalue p.d.f. is proportional to (1^), and the corresponding n- 
point distribution function is given in terms of the orthonormal polynomials {pj(a:)}j=o,i,2,... 
associated with the weight function W2{x) according to the formula 

Pn{xi, . . . ,Xn) = det[KN{Xi,Xj)]-^^^..^^ , (2.3) 

where 

N-l 

KN{x,y) = [w2{x)w2{y)]^^'^ ^Pi{x)pi{y) • (2.4) 

1=0 

Substituting into (p.2|), the Fredholm theory of integral operators then gives 

E{Q;I) =det(I-KAr) , (2.5) 

where IKat is the integral operator with kernel K]\j{x,y) defined on the interval /. A crucial 
point is that (2.4) can be summed according to the Christoffel-Darboux formula and so 
written in the special form [12| 

x-y 

where with denoting the coefficient of x^ in pn{x) 

(j){x) = ( — —W2{x) ) Pn{x) , 

\ aN J ^ ^ 

(2.7) 

f a-N 1 \^ 

Ipix) = ( —W2{x) ] PN~l{x) ■ 

\ aN J 

We suppose furthermore that (j){x) and 'ijj{x) satisfy the recurrence-differential relations 

ra{x)(j)' [x) = A[x)(f)[x) + B{x)ip{x) , 

(2.8) 

m{x)ip (x) = —C{x)(j){x) — A{x)^{x) , 

where the coefficient functions m{x), A{x)^ B{x),C{x) are polynomials in x. 

In the above setting Tracy and Widom |jl^ have derived a set of coupled differential 
equations for the determinant ( |2.5D , as well as some auxiliary quantities. The latter we 
introduce in the following definitions. 
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Definition 1 Let A = A{x,y) denote that the integral operator A has kernel A{x,y). Then 
the kernels p{x, y) and R{x, y) are specified by 

{l-K)-'=p{x,y) , 

(2.9) 

K{l-K)-' = R{x,y) , 
The operator K{1 — K)^^ is called the resolvent and R{x,y) the resolvent kernel. 
Definition 2 For k € Z>o the functions Qk and Pj. are defined by 

Qk{x) = j dy p{x,y)y^4>{y) , 

' (2-10) 
Pkix) = j^dy p{x,y)y'''ip{y) , 

and their values at the endpoints aj of I are denoted qkj,Pkj so that 

Qkj = Qk{aj) = lim Qk{x) 

x—^aj 

Pkj = Pk{aj) = lim Pfc(x) 

Where there is no confusion, we denote qoj,Poj by qj,Pj- 
Definition 3 The inner products u, v, w are defined by 

u = (</)|Q) = [ dy Qo{y)m , 



(2.11) 



V = {i;\Q) = j^dy Qoivmy) = (0|P) = j^dy PoiyMy) , (2.12) 

w = {ij\P) = j dy PQ{y)^{y) . 

Our goal is to characterise the probabilities (^]J) and ( |1.5D as the solution of certain non- 
linear differential equations. Following |18| this is achieved by specifying partial differential 
equations for the quantities qj,pj,u,v,w and R{aj,ak)- These equations come in two types: 
a set of universal equations which are independent of the recurrence-differential equations 
( |2.8D , and a second set of equations which depend on the details of ( |2.8D . Let us first present 
the former. 

Proposition 1 For general functions (l){x),ip{x) we have the relations 

^logdetfl -i^) = (-iy-^R(ai,ai) , (2.13) 



and for j ^ k, 



and 



N qjPk-PjQk 
R{aj,ak) = — , (2.14) 



d 

- — R{aj,aj) = {-l)^R{aj,ak)Riak,aj) , (2.15) 
oak 
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with 



dqj_ 

dak 

dpj_ 
dak 



{-lfR{aj,ak)qk , 
{-lfR{aj,ak)vk , 



(2.16) 



and 



du 
dak 



(-1) 



k„2 



— = {-ifqkPk , 
oak 

— -( 

dak - ^"'^ • 



(2.17) 



The second set of equations, which depend on the details of ( |2.8| ), give the j = k cases 
of (|1|) and ( ^l6|) . Now for weights (|1.2| ) the equations hold for m(x) a quadratic and 
A{x) , B (x) , C (x) linear functions, and thus of the general form 



m(x) = /io + fJ-ix + ii2X^ , 
A{x) = ao + aix , 
B{x) = /?o + , 
C{x) = 70 + 71 a; • 



(2.18) 



One then has the following equations [18|. 



Proposition 2 In the case that (l){x),ip{x) satisfy the equations ( ^-ai) with coefficient func- 

dqi 



tions ( 2. It ) we have 



nii 



dai 



rui 



dai 



[ao + aiai + 7iii - I3iw - H2v]qi 
+ [Po + Piai + 2aiu + 2Piv + fi2u]pi 

2M 

-'^{-l)''R{ai,ak)qkmk , 

[-70 - 71 Oi + 271 -u + 2aiw - fi2w]qi 
+ [-ao - aittj + Piw - 71 li + ^i2v]Pi 

2M 

-'^{-l)''R{ai,ak)pk'mk , 

k^i 



(2.19) 
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miR{ai, ttj) = [70 + 7iai - 2711; - 2aiw + fi2w]qf 

+ [Po + PiUi + 2a\u + 2(5iv + ^Ji2u]p'l 

+ [ao + aitti + 7in - - fj,2v]2qiPi 

and furthermore, 

d 

— [miR{ai,ai)\ = 2aiqiPi + PiPi + JiQi 

2M 

-^{-l)''mkR\ai,ak) , (2.21) 



where mi = m{ai). 
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3 The Gaussian Ensemble 



In this section we take up the problem of computing ( |1.5D . We recah (see e.g. |T^]) that for 
the Gaussian weight the corresponding orthonormal polynomials {pn{x)}n=o,i,... are 
given in terms of the standard Hermite polynomials H]\f{x) by 

Pn{x) = 2-^/2vr-V4(^r!)-i/2^^(^) ^ (3.I) 
and for the coefficient a at of in pn{x) we have 

Substituting in ( p.7lJ2.8| ) and making use of the differential relation for the Hermite polyno- 
mials gives Ijl^ 

(j)'[x) = -x(t>{x) + V2Nip{x) , 

(3.3) 

^'(x) = - V2N(p{x) , 

and so we have 

m{x) = 1, ai = -1, /3o = 70 = \/2iV, ao = /3i = 71 = . (3.4) 
3.1 Differential Equation 

The interval from which the eigenvalues are excluded in the probability ( |1.4| ) consists of 
the union of the two intervals (— 00, — s) and (s,oo). Thus in ( |2.lD we have M = 2, and 
ai = —00, 02 = — s, as = s and 04 = 00. Because ai and 04 are infinite, only the quantities 
in Propositions 1 and 2 relating to 02 and 03 are of interest. Furthermore it turns out that the 
evenness symmetry of / and the odd/even symmetry of ( ^.iD implies the equations relating to 
a2 are equivalent to those relating to 03. Thus <j){—x) = {—\)^(f){x), Tp{—x) = {—l)^~^'ip{x), 
and so KN{—x,—y) = KN{x,y), which in turn implies p{—x,—y) = p{x,y) and so q2 = 
(-1)^^93, P2 = R{-x, -y) = R{x, y). 

Proposition 3 The coupled differential equations for q = q^^p = p^ and R{s),R{—s,s) of 



the finite N GUE on the interval I = (—00, — s) U (s,oo) are 

^lnE2 = 2R{s), (3.5) 
as 

q' = -sq + p[(3o - 2u] + 2q^p/s , (3.6) 

p' = +sp - q[jo + 2w] - 2qp^/s , (3.7) 

u = -2q^ , (3.8) 

w' = -2p^ , (3.9) 

R{-s,s) = {-l)^-'qp/s , (3.10) 

R{s) = q^bo + 2w] +p2[/?o - 2u] - 2sqp + 2gV/s , (3-11) 

R'{s) = -2qp - 2^V/s^ , (3.12) 
where primes indicate derivatives with respect to s. 
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Proof - The first equation follows from ( 2.13|), (|3 .6|) and 



,, D) follo'^ fro m (|2.19| ), ( |3.8|) an d 

(|]|) follow from 1^1% , (|3.10D follows from (|2.14D , (^.llD from and (|3.12|) from (|2.2lD . 

In the derivation of ( |3.5| ) and ( 3.12| ) use is made of the general formula for the total derivative 



d 
ds 



fi{s,s) 



_d_ 

da? 



f I {as, as) 



_d_ 
da-. 



fiia3,a3] 



03 = — 12=S 



(3.13) 



valid whenever fi{a2,a2) = //(as, 03). □ 

There are a number of differences when one compares this set of coupled ODEs with that 



arising from a single interval (s, 00), as described in Ref.[18|. Firstly there is the appearance of 
factors of two, reflecting the contributions from the two free endpoints, but more significantly 
is the addition of nonlinear quadratic terms in |3.7|, |3.11 , |3.12|). Regarding the reduction 
of these equations, it is straightforward to show that ( |3.6D and (l3?7|) together imply 



or equivalently 



pq = uw- y2/3ou^ + y27o^i , 



Polo - ^pq = iPo - 2n)(7o + 2w) 



(3.14) 



(3.15) 



which is useful in the subsequent determination of a single equation for R{s). In fact it is 
possible to reduce the coupled system of PDEs in Proposition 3 down to a single ODE for 
R{s), or a single ODE for the quantity 



R{s) = (-l)^i?(-s,s) . 



(3.16) 



The boundary conditions satisfied by R{s) and R{s) are found from the large s form of 
the kernel 



R{s) 



S—fOO 



(3.17) 



which implies 



Ris) 



R{s) 



1/22-N 

[H%{s) - Hr,+,{s)Hr,.,{s)] 



7r-V22- 



(3.18) 



(iV-1)! 



Of course we could replace the Hermite polynomials by their leading order term, however of 
subsequent interest will be a scaling limit which requires the full expression as written. 

We can reduce this coupled system of PDEs into a single second order ODE, and give 
two possible forms for this. 

Proposition 4 The system of ODEs in Proposition 3 is equivalent to the following second- 
order differential equation for R{s) 



sR" + 2R' + 8NsR + 2As'^R^ 



A[s + 2R]'' {R + sR'f + 8Ns'R' + l6s^R^ 



(3.19) 
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or to the second-order differential equation for R{s) 



sR" + 2R' = 2s[s -h]- 2h^J{R + sR'Y - 4s2[s - h]R - 2Ns'^[s - hf 
where we denote 

h = Vs^ - 2R' . 



Proof - We adopt the notation a = sR, b = sR for simplicity. Using ( |3.6[ ) and (| 
with the definition ( 3.16| ) we find the relation 



b' = q\jo + 2w)-p\(3o-2u) , 
and utilising ( 3.11| ) and ( |3.12D the corresponding relation for a is 

a' - Asb = q^{-fo + 2w) +p^[Pq - 2u) . 



(3.20) 

(3.21) 
together 

(3.22) 
(3.23) 



By squaring and subtracting the right hand sides of these two equations and employing the 
integral ( p. 15 ) to eliminate the cross term we have 



{a'f = Sab + 8Nb^ + {b'f . 
There is also another relation which follows directly from ( |3.12| ) and reads 



sa 



' - a + 2s'^b - 26^ . 



(3.24) 



(3.25) 



If one proceeds to eliminate a using ( ^.24 ) and ( 3.25 ) then a second-order differential equation 



is obtained for R, namely (3.19). However if one eliminates b then a more useful second order 
equation for R emerges 



1/2 [s[s T Vs^ - 2R'] - i/2a"}' = 

(s^ - 2R') {1/2(0')^ - 2s[s T Vs^ - '^R']a - Ns^[s T \/s2 - 2i?']2} . 
Firstly the square-root appearing above is well defined as 

- 2R! = {s- 2Rf . 



(3.26) 



(3.27) 



One now faces a decision regarding the choice of the branch to be taken for the square root 
of the right-hand side of ( 3.26| ), and if we take the negative branch then the correct choice is 
given in ( 3.2C| ). This choice is the only acceptable one given the boundary conditions ( p. 18 ) 
which impose that all quantities a, R, R and derivatives vanish exponentially fast as s ^ 00. 
If one takes the upper, negative sign then s — h also vanishes in the same way and is consistent 
with the ODE ( 3.26| ) but if the other choice is made then s + h = 0(s) in this limit and 
is incompatible with the ODE, both in the leading order s-dependence and the sign of the 
coefficient. □ 

3.2 Special Cases 



Here we calculate (|l.4| ) in the cases = 1 and N = 2 directly from the expansion (|2.2| , 2.3 , 



|2.6[ ), which truncates after term n = N . The formulas ( |3.5[ ) and (3.27) then allow R{s) and 
R{s) to be computed in terms of the error function. 
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Proposition 5 Consider the probability (I.4) and the associated functions R{s) as specified 
by (3^) and R{s) as specified by ( 3.16 ). For N = 1, with the the standard definition of the 
error function erf(s), we have 

^2(0;/) = erf (s) , 



R{s) 



7r^/2erf(s) ' 



(3.28) 



R{s) 



while for N = 2 



7r^/2erf(s) ' 

£2(0; I) = erf(s) ferf(s) - 271-^^"^ se'' 
Ris] 



+ 



2 

s^e 



7rV2erf(s) i/aTrVSerf (s) 



se 



(3.29) 



Ris) 



2 -s^ 



7rV2erf(s) i/avrVaerf (s) 



se 



Using the kernel above, and ( |2.7D and (^2), the results follow from direct evaluation of the 
integrals. □ 

We can check (by hand and using computer algebra) that these expressions satisfy the 
appropriate equations ( p. 20 , 3.191) in Proposition 4, and furthermore satisfy the appropriate 
boundary conditions s ^ 00, equation ( ^.18 ). 

For small s we find from ( 3.28 , 3.29]) that 



2s 3 45 



16 



945 



14175 



-s' + 



32 



-s^ + 



1472 



-,11 



93555 638512875 



R 



N=2 



2 14 
--—s + 



248 



128 



51104 



s 15 1575 



23625 27286875 



s' + 



1356032 



„9 



6898816 



41 



5320940625 558698765625 



(3.30) 



We can also determine the small-s expansion for general A^. The form of this expansion can 
be understood from the small s behaviour of (|1.4D. Now 



E2{0; (-00, -s) U (s, 00); e-^' ; ^) = 77 H / ^^^e"^? J] '^'^ " ^. 

1=1-^-^ j<k 



1=1 

N2 N 1 



1=1-^ j<k 



~ aos +028 ^ + ... , 



(3.31) 

where the final equality follows by expanding the exponential. According to ( ^.51) and ( [3.27]) 
this implies 



R{s) ~ — + 5is + 635^ + . . . . 
s 



(3.32) 



Making an ansatz of this form we find that all the coefficients are uniquely determined by 
( 3.20| ). Use of computer algebra gives the expansion in the following result. 
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Proposition 6 For fixed N, the asymptotic expansion of R{s) for small s is given by 

N{2N^-1) 2N'^{4:N^ -9N^ + 3) ^ 
^""^ 4N^-1 ^ (4A^2_ 1)2(4^2 _ 9)^ 

8iV3(47v4 _ i3Ar2 + 6) g 
^ (4iV2 _ i)3(4Ar2 _ 9)(4Ar2 _ 25) * 

8A^2(i28A^^° - 1312iVS + 3304A^6 - 3430iV^ + _ 315) ^ 

^ (4Af2 _ i)4(4Ar2 _ 9)2(4Ar2 _ 25)(47V2 - 49) ^ ' ^^'^^^ 

It is of interest to note that in fact the only Laurent series about s = which satisfies 
( |3.20| ) has the structure ( 3.33 ). To see this it is simplest to consider the equation ( 3.19 ). 
Making the Ansatz 



00 



sR = Y^ rns"- , (3.34) 



n=k 



we find that with a lower exponent k and k > 1 the only solution for the coefficients is the 
null solution. However if k < —1 then one finds = . . . = r_i = 0, and furthermore the 
relations obtained by equating the coefficients of s~'^, s^^, . . . are the same as for the case of 
k = 0. Explicitly equating coefficients of s"^ in the ODE yields 

rl + 8Nrl + 16r|] = . (3.35) 

Now ro 7^ 0, therefore the odd-index terms must vanish because sR{s) is even, as is seen 
from the ODE ( p. 19 ). Thus ri = and an unique solution is found for ro = —N/2. The 
coefficient of automatically vanishes while equating coefficients of gives 

iV2 

Choosing the former root, we find by equating coefficients of higher powers of s that each of 
rs, r4, . . . is now uniquely determined. 



3.3 Edge Scaling 

To leading order the support of the GUE is within the interval {—V2N, V2N). The statistical 
properties of the eigenvalues in the vicinity of the edge (A ~ \/2N say) can be studied by 
introducing the scale 

A ^ V2N + , , (3.37) 

which for N ^ 00 makes the separation between the largest and the second largest eigenvalue 
of order unity. In ( |l.4| ) the probability that the region —00 to the vicinity of the lowest 
eigenvalue, and the highest eigenvalue to 00, is free of eigenvalues can be studied in this limit 
by writing 

s = V2N + , . (3.38) 
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Because the two intervals in become independent in this hmit we would expect 



^lim E,iO; i-oo, -(Vm + -^)) U (Vm + oo); ; N) 



^/2iVi/6 ' 



[Ef\0;{t,^))y , (3.39) 



where S|°^*(0; (t, cxd)) refers to the probability that the edge (t, oo) is free of eigenvalue in the 
infinite GUE scaled according to ( |3.37 ). Now it follows from (^]^) that 



/ POO 

E2iO;{-oo,-s)U {s,oo);e'^^;N) = exp (-2 / duR{u] 
On the other hand it is known that [ p^ ] 

Ef\0; (t, oo)) = exp duR'°^\u) 

where W°^^{t) satisfies 



iSoft 



R 



soft 



^j^soit _j_ ^soft 



, 



(3.40) 



(3.41) 



(3.42) 



which is the Jimbo-Miwa-Okamoto form of a particular Painleve-II equation, subject to the 
boundary conditions 



R^oit^^-j ^ _t[Ai(t)]2+ [Ai'(t)l^ . 

t^oo 



Thus for ( 3.39| ) to be valid we must have 

1 



lim 



:R(V2N + 



) = R'"^\t) 



The equation (|3.42| ) can be verified by putting 



r{t) 



^ R{V2N+ * 



(3.43) 



(3.44) 



(3.45) 



and showing that for N ^ oo r{t) satisfies the differential equation ( 3.42| ) and the boundary 
condition ( p.43| ). Regarding the latter point, note that ( |3.17 ) implies 



rit) 



(3.46) 



where Kn{x, y) is specified by substituting (|3.lD and (|2.7|) into (|2.6|) . But we have the known 
limiting behaviour |10| 



lim ^ 



/2N + 



Y 



1/6^ 



Jv — Y 



which immediately imphes ( [3.43| ). Thus it remains to show that for N ^ oo ( [3.45| ) satisfies 

(H). 
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In fact it is more convenient to work with the ODE (3.19), and introduce the scaled 
quantity 



V2 



(3.48) 



(this choice of scaHng is consistent with the relation ( |3.12D between R{s) and R{s)). With 
this substitution one finds that for ^ cxd ( p.l9| ) reduces to 



2ff - (f=)^ + 16f^ - 4tf^ = , 
while the relation ( 3.12| ) reduces to 

f = 2r . 

Substituting ( p. 501 ) in ( p.49| ) gives the third order ODE 



.=.2 







(3.49) 



(3.50) 



(3.51) 



Differentiating (3.42) and using the original equation once more, it is easy to reduce it to 
( 3.51 ), thus establishing the required result. Alternatively one prove the same correspondence 
by making the substitution 



(3.52) 



in the ODE (3.49), and then show this leads to a Painleve-II equation with a = 0, namely 

q = 2q^ + tq . (3.53) 
This equation has been derived for the kernel ( p. 47 ) in [p^). 
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4 The Jacobi Ensemble 



In this part of our study we will treat the Jacobi ensemble, for general a, /3 in some parts, 
but mostly we will consider only the symmetrical case a = (5 with our objective to compute 
the quantities (|l.5| ). For the Jacobi weight in ( |l.2| ) the orthonormal polynomials are given in 
terms of the standard Jacobi polynomials by 



PNix) 



2N+a+P+l N\r{N+a+P+l) 

2*^1 r{N+a+i)r{N+p+i) 

with the corresponding coefficient of such that 

N{N+a)iN+l3)iN+a+l3) 



1/2 



P, 



N 



(x) 



QjV-1 



(2A^+a+/3)2(2iV+a+/3+l)(2iV+a+/3-l) 



1/2 



(4.1) 



(4.2) 



Making use of the differentiation formula 



{2N+a+P){l-x^)-^Pi^'''\x) = N[a-P-{2N+a+p)x]P^^'''\x) 
ax 

+ 2{N+a){N+P)P^^l^,\x) , (4.3) 



and the three term recurrence for the Jacobi polynomials gives that the coupled first order 
equations hold with |T8[ 



m{x) = 1 — , 

, /32-a2 2N+a+f3 

Alx) = , —- X = an + OL\x 

^ ' 2(2A^+a+/?) 2 



^, . 2^Af(Af+a)(Af+/3)(Af+a+/3) /2A^+a+/?+l ^ (4.4) 
^^""^ 2iV+a+/3 V2iV+«+/5-l " ^° ' 



^ 2v^A^(A^+a)(A^+/3)(A^+«+;g) / 2A^+a+/3-l _ 



4.1 Differential Equations for End Intervals 

We will begin by considering the first probability in ( |1.5[ ), generalised so that the weight 
function is the general (nonsymmetric) Jacobi form in ( |1.2D . We shall adopt the conventions 
93, P3 = Q2,P2 = and R+ = R{s,s),R- = R{-s,-s),Ro = R{-s,s). 

Proposition 7 The coupled differential equations for the finite N JUE on the interval 
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(— 1, — s) U {s, 1) for general a, (3 are 
[\nE2]' = R-+R+ , 



Rq 



q+p- - Q-P+ 

2s 

[70 - w{2ai + l)]ql + [po + u{2ai-l)]p'i + [qq - ais + v]2q_p. 
+ 2s{l-s^)Rl , 

[7o - w{2ai + l)]ql + [Po + n(2ai-l)]p+ + [uq + ais + v]2q+p^ 
+ 2s{l-s^)Rl , 



V = -q^p^ - q+p+ 



w 



2 2 
-P- - P+ , 



{l-s')q 



-[ao - ais + - [/?o + u{2ai-l)]p- - 2(l-s^)g+i?o , 
-[-70 + w{2ai + l)]q- - [-ao + "i^ - v\p^ - 2{1-s^)p+Rq 



{l-s')q'+ = [ao + ais + + [^o + u{2ai-l)]p+ - 2(l-s^)g_i?o , 



(1-s )p+ = [-70 + w(2ai + l)]g+ + [-ao - ais - v]p+ - 2(l-s )p_-Ro , 
[{l-s^)Rj\' = -2aig_p_ - 2{1-s')RI , 
[{l-s^)R+]' = +2aiq+p+ - 2{1-s^)rI . 



(4.5) 
(4.6) 

(4.7) 

(4.8) 
(4.9) 
(4.10) 
(4.11) 
(4.12) 
(4.13) 
(4.14) 
(4.15) 
(4.16) 
(4.17) 



The first equation follows fro m (|2.13[ ) , ( |2.5| ) and the next equality in ( [4.6[) follows from ( |2.14D , 
(13) and (U) follow from (|2.2C| ), while (|4.16D and (p^ ) follow from (|2.2lD . Furthermore 
(I^-(pID follow from ( |2T7D and the first equality in (j^ and (p^)-(^l5|) follow from 
§A^ ) and (0). □ 

We note from ( p^ , (|4T^ ) and ( pO|) the integral 

(l-s2)(i?_ - i?+) = 2aiz; . 
The boundary conditions satisfied by R{s, s) as s — > 1~ or s - 
Nir{N+a+p+l) (l-s)"-i(l+s)^-i 



(4.18) 

-1+ can be expressed as 



R{s,s) 



r{N+a)r{N+P) 2«+/3(27V+a+/3) 

P^-a^ 

(2iV+a+/3)(2iV+a+/3+2) 



2{N+l){N+a+P+l) ia,P) ia,(3). . 



+ 



2N+a+P+2 

2N{N+a+P) 
2N+a+p 



N 



(4.19) 



Our objective is to use equations such as ( 4.18 ) to reduce the equations of Proposition 
7 down to a single equation for R-^-{s). For this purpose we restrict attention to the case 
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a = (3, then we see from (|4.4| ) that oq = 0, while the fact that 

4"'"^(-^) = (-l)M"'"^W> (4-20) 

imphes f = 0, ^2 = (— l)^93i P2 = (— 1)^~^P3 and i?+ = i?_ (= i? say), thus reducing the 
number of unknowns in Proposition 7. 

Proposition 8 The coupled differential equations for the finite N JUE on the interval 
(—1, — s) U (s, 1) for P = a are 

[lnE2]' = 2R , (4.21) 

R{-s,s) = (-l)^-i^ ^ (-l)^-^i?o , (4.22) 

(l-s2)i? = - u;(2ai + l)]g2 + [/Jq + u{2ai-l)]p'^ + 2aisgp + 2s(l-s2)i?2 ^ (4.23) 

= -2g2 , (4.24) 

w' = -2p^ , (4.25) 

2Cl_ q2N 

(l-s2)g' = +aisg+[/3o + ^,(2ai-l)]p + ^ Lq'^p , (4.26) 

s 

2('l_s2\ 

{l-s'^)p' = -aisp + [--fo + w{2ai + l)]q - ^ ^-qp"^ , (4.27) 

s 

[{l-s')R]' = 2a,qp - ^^^^^qV , (4.28) 

where we have redefined Rq from the previous usage (Rq i— > (— l)^~^i?oJ and made the 
notation q+ = q,p+ = p. 

We now indicate how to reduce such a system to a single third order differential equation 
for R = R{s) = R{s,s). 

Proposition 9 The coupled set of ODEs given in Proposition 8 reduce to the third order 
ODEfora{s) = {l-s'^)R{s), 

il-s'fa'" - 2s{l-s^)a" - 2s-\l-s^)a' - {l-s^f{a"f /a' 

S aiS + H ,, 1 n 2 

+ "^2^^^[(^-'^^ ^^"^ ^] 

2 

+ AaiHa + 2[ais + H] (ais^^ - Ha') + 2-^H\H - aisl = , (4.29) 



where H = \J a\s'^ — 2(1 — s^)^'. 

Proof - Firstly we make some auxiliary definitions similar to the Gaussian case 

a = qp , 

^ = [/3o + n(2ai-l)]/ , (4.30) 
B = [-7o + u;(2ai + l)]g^ . 
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Combining the relations ( [4.26 ) and ( |4.27| ) we find 

{l-s^)a' = A + B 



(4.31) 



and we differentiate this once more, using the relations iWM) , (HI) (|436|) and (|07|) , to 
arrive at 



{l-s^fa" =2s{A + B)- Sai{l-s^)a^ 

+ [2ais + As~^{l-s^)a]{B - A) + AAB/a 



(4.32) 



Now the idea is to express A, B in terms of a and its derivatives, so we employ ( 4. SI] ) above 
and the relation 



sa' + a — Aaisa = A — B 



(4.33) 



which fohows from ( [4. 231) and ( [4.28| ). In order to express a in terms of o"'(s) we have to solve 
the quadratic relation ( |4.28D for a and this is how the square-root variable H arises. This 
quantity is well defined because 



a\s^ - 2(l-s^)(T 



2\ / 



2 a — ais 



(4.34) 



Having expressed A, B and a in terms of a and it derivatives it is then a matter of substituting 
these into ( |4.32| ) and after considerable simplification we find the final result ( }4.29| ). □ 



4.2 Special Cases for End Intervals 



In this part we present the calculations for the first two finite- A'^ cases, that is iV = 1 and 
= 2, by direct means using the probability E2{0;I). 

Proposition 10 The probability that there are no eigenvalues in the interval I = (—1, — s)U 
{s, 1) of the JUE with N = 1, and general a ^ [3 is 

E2{Q;I) = /(i+.)/2(a+l,/5+l) - /(i-s)/2(«+l, /3+1) (4.35) 

and for N = 2, is 

^2(0;/) = (a+/3+3) [/(i+.)/2(a+l, /3+2) - /(i_,)/2(a+l, /3+2)] 

X [/(i+,)/2(a+2,/3+l) -/(i_,)/2(a+2,/3+l)] 
- (a+/3+2) [/(,^^)/2(a+2,/3+2) -/(i_,)/2(a+2,/3+2)] 

X [/(,+,)/2(a+l,/3+l)-/(i-s)/2(a+l,/3+l)] (4.36) 
where the normalised incomplete beta functions are defined by 

in terms of the incomplete and complete beta functions, Bx{a,b) and B{a,b) respectively. 
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Proof - as in the case of Proposition 5, this follows from the expansion (|2.2[) with (|2.3[) 



which truncates after the term n = N used in conjunction with (|2.6D , (2.7), ( [4.1| ) and ( |4.2D 
Extensive use of the identities for the incomplete beta function 

Ix{a,b) = 1 - a) , 

{a+b)Ix{a, b) = alx{a+l, b) + blx{a, b+\) , 

are made to reduce the number of their occurrences. □ 

In the symmetric case a = (3 considerable simplification ensues and we have the following 
results. 

Proposition 11 In the case (3 = a, the probability ^-35 ) and the associated quantities of 
Proposition 10 are given in terms of the Gauss hypergeometric function 2-^1(0, 6; c; z) by 



^2(0;/) 



ais) 



^ (l-g2)a+l 

'/'s2Fi(-a,i/2;3/2;s2) 



S2Fi(-a,y2;3/2;s' 



(4.38) 



H{s) = {a+l)s + 



2\a+l 



and similarly for the probability ( j.St ), 
2a+3 



E2{0;I) 



ais) 



42"+iS2(a+l,a+2) 



2Fi(-«, 1/2; 3/2; s') 2Fii-a, 3/2; 5/2; s' 



(l-s^) 



2\a+l 



2s 



1 



+ 



2Fi(-a,i/2;3/2;s2) 2Fi{-a,Y2;V2;s^)j ' 



(4.39) 



H{s) = {a+2)s + 



2\a+l 



1 



+ 



2Fi{-a, 1/2; 3/2; s2) ^Fii-a, 3/2; SA; ^2 



Proof - These follow from the reduction of the incomplete beta functions in the symmetric 
case to Hypergeometric functions, such as 

^(i+s)/2(a+l,a+l) -/(i_,)/2(a+l,a+l) = 4^;g(^^f^^^ 2i^i(-a, 1/2; 2/2; s^) , (4-40) 

the differentiation formulae for these Hypergeometric functions, like 

^[s2Fi(-a,i/2;3/2;.2)] = (i_s2)" (4.4I) 

and the use of their contiguous relations 
(2a+3)i/3s32Fi(-a,3/2;5/2;s2)_5 2Fi(-a,i/2;3/2;s2) + s2Fi(-a-l,3/2;3/2;s2) = . (4.42) 

□ 

One can show that these two specific cases are solutions to our third order differential 
equation ( [4.29 ), after noting the contiguous relation above ( [4.42 ) linking the two Hypergeo- 
metric functions in the case of = 2. 
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Another special case for which the probabihty can be computed independently of the 
ODE ( [4.29|) is that with q = /? = and general N. We then have 



E2{0;{-h-s)U{s,iy,X[o,i];N) = ^ r dXi--- f cIXn J] I^j " ^^1' ' (4-43) 

-^'^ l<j<k<N 

where C is such that £'2 = 1 for s = 1, which can be evaluated by a change of variable. 

Proposition 12 The probability o,nd the associated quantities of Proposition 10 have 

the evaluation 



'yis) = N'^^^ > (4.44) 
H(s) = - , 



n 



This exact form of o"(s) can be verified to satisfy ( 4.29| ) with a = /? = 0. Furthermore 



we can take the limit q ^ in ( [4.38 ) and ( 4.39 ) and reclaim ( [4.44 ) in the cases = 1 and 
N = 2 respectively. 

4.3 Jacobi to Hermite Limit for the End Intervals 

From the definition 

E2{Q■,{-l,-s)u{s^)■,{l-\^r■N) = 

N 



^ r dX,... r dX^fl{l-Xf)- II \Xk-XjW (4.45) 



1=1 l<3<k<N 



where again the constant C ensures the normalisation. Replacing s by t/^/a^ changing 
variables Xi^ Xi/ ^/a and taking a ^ 00 shows 



£;2(0; (-1, -tj^) U (t/V^, 1); (l-A^)"; iV) ~ 



a— »oo 



E2{^\ (-00, -t) U (t, 00); e"-^ ; N) . (4.46) 
Recalling ([3.5[ ) and ( [4.21 ) this is equivalent to the statement that 



lim ^i?(^) = i?GUE(^) ^ (4 47) 

where R{s) on the left-hand side is as in Propositions 7 and 8 while R^^^{t) is the R defined 
in Proposition 3. Thus it must be that with 

1 i?(^) = rit) , (4.48) 



in ( 4.2S| ), taking the limit a — > cxo must take the third order equation into an equation 



equivalent to the second order equation ([3.2C| ) 
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To verify this, we first note that the substitution (4.48), to leading order in c? reads 



r'" + 2 



+ 2r'h[t - h] 



2 2 

2 



(4.49) 



and that differentiation of the second order Hermite ODE, ( 3.20| ), gives 
tr'" + 3r" =4t + Ah^{tr' + r + 2Nt) - 2h 

- 8t{r + Nt)h -{t- r"){2t^ - tr" - 2r')/h^ 



(4.50) 



with h = Vt^ — 2r', as before. By effecting a suitable subtraction of these two equation in 
order to eliminate the term in r'", and using the second order ODE once more, one can show 
this difference is identically zero. 

Employing this limit one can show that the special cases of A'^ = 1 and N = 2 for the 
Jacobi weight ( 4.3§| ) and ( 4.39 ) lead exactly to those for the Hermite weig hts (I^H) and 
( |3.29| ), respectively. 



4.4 Differential Equations for an Interior Interval 

We will now consider our last case, the second probability in ( |1.5| ), in a parallel manner 
to that of the previous case. Treating first the nonsymmetrical form a ^ (3, we make the 
new conventions q2,P2 = <1+,P+, Qi,Pi = <1-,P- and = R{s,s),R- = R{—s,—s),Ro = 
R{-s,s). 

Proposition 13 The coupled differential equations for the finite N JUE on the interval 
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(— s, s) for general a, (3 are 



{l-s^)R- 

{l-s'^)R+ 

u' 
v' 
w' 



— -R_ — Rjf- , 

q+p- - Q-P+ 

2s 



[70 - w{2ai + l)]ql + [f3o + u{2ai-l)]p'i + [oq - ais + v]2q^p^ 



2s{l-s^)Rl , 



[70 - w{2ai + l)]q1 + [f3o + u{2ai-l)]p1 + [oq + ais + ^;]2g+p_, 



+q- + ql , 



+q-p- + , 

+p- + pi , 



-[qo - ais + - [/?o + u{2ai-l)]p^ + 2(l-s^)g+i?o , 
'[-70 + w{2ai + l)]q- - [-ao + otis - v\p- + 2{l- s^)p+Rq 



{l-s^)q\ = [ao + ais + + [/3o + n(2ai + 2(l-s^)g_i?o , 
(l-s^y^ = [-70 + u;(2ai + l)]g+ + [-qq - qis - v]p+ + 2{l-s^)p-Ro , 
[{l-s^)R-]' = -2aig_p_ + 2{1-s'^)RI , 
[il-s^)R+]' = +2aiq+p+ + 2{1-s^)RI . 



(4.51) 
(4.52) 

(4.53) 

(4.54) 
(4.55) 
(4.56) 
(4.57) 
(4.58) 
(4.59) 
(4.60) 
(4.61) 
(4.62) 
(4.63) 



Proof - these follow in an entirely parallel manner as for the derivation of the previous set 

(iD-(ii3). □ 

The only, apparently minor, difference between the interior interval set of equations, 
(^^)-(|46^), and the endpoint interval set, (|45|)-(p^), is a change in sign of a number of 
terms in the expressions for derivatives. However, in reality, the two cases are quite distinct. 

Again we note that ( [4.62 ), ( |4.63| ) and ( [4.56| ) imply the integral 



[l-s'){R, 



R- 



2aiv . 



(4.64) 



The boundary conditions satisfied by R{s, s) now apply as s ^ and the limiting value takes 
the same form as ( 4.19D . 

We continue by considering the symmetrical case a = /3, and find again that oq = 0, 
and the parity relation ( 4.20| ) implies v = and R-^- = which is denoted by R. 



Proposition 14 The coupled differential equations for the finite N JUE on the interval 
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(— s, s) for P = a are 
[\nE2]' = -2R , 

i2(-.,.) = (-l)^-i^ = (-l)^-% 
s 



u' = +2g2 , 



(l-s^)g' = +aisq + [Po + u{2ai-l)]p 



{1-s )p = -aisp + [-^Q + w{2ai + l)]q + 
[il-s')R]' = 2a,qp+'^-^^^^qV . 



2{l-s^) 



-qp 



(4.65) 
(4.66) 



{l-s^)R = [70 - w;(2ai + l)]g^ + [/?o + u{2ai-l)y + 2aisqp - 2s{l-s')Ri , (4.67) 



(4.68) 
(4.69) 

(4.70) 
(4.71) 
(4.72) 



where we employ the notation q2 = q,P2 = P- 

Again such a system can be reduced to a single third order differential equation for R{s). 



Proposition 15 The coupled set of ODEs given in Proposition 13 are equivalent to the third 
order ODE for a{s) = {l-s'^)R{s), 

{1-s^fa'" - 2s{l-s^)a" - 2s-\l-s^)a' - {l-s^f{a"f /a' 
s ais + G r 2x // o -1 /i2 

2 

+ 4aiGa - 2[ais + GUais''^ + Ga') + 2-^G[G - ais] = , 



(4.73) 



where G = afs"^ + 2(1 — s^)^' 



Proof - this proceeds in an identical manner to the proof of Proposition 9 but with a number 
of minor alterations. As in the previous reduction we have to express the following equation 



{l-s^fa" =2s{A + B) + 8ai(l-s2)a2 

+ [2ais - 4:S~\l-s^)a]{B - A) + AAB/a 



(4.74) 



in terms of a and its derivatives alone. We have an analogous expression for a' which is 
quadratic in a and its solution contains the square-root variable G defined above. This is 
well defined due to the relation 



a{s^ + 2{l-s'^)a' 



2 a + ais 



(4.75) 



The final result is the above third order ODE ( 4.73 ). □ 
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4.5 Special Cases for the Interior Interval 

In this part we present the analogous results for the first two finite- cases = 1 and N = 2 
by direct calculation of the probability £"2(0;/). 

Proposition 16 The probability that no eigenvalues are found in the interval {—s, s) for the 
JUE with N = 1, and general a ^ (3 is 

^2(0;/) = l-/(i+,)/2(a+l, /?+!)+ /(i_.)/2(a+l,/3+l) (4.76) 

and for N = 2, is 

£2(0;/) = (a+/3+3) [l + /(i_,)/2(a+l,/3+2)-/(i+,)/2(a+l,/3+2)] 

X [1 + V«)/2(«+2, - ^(i+s)/2(«+2, 

- {a+P+2) [1 + /(i_,)/2(a+2, /3+2) - /(i+,)/2(a+2, /3+2)] 

X [i + /(i_^)/2(a+l,/3+l)-/(i+,)/2(a+l,/3+l)] (4.77) 
in terms of the normalised incomplete beta functions, Ix{a,b). 

Proof - these follow from direct evaluations of the terminating series for the probability 
E2{0; I) given in (2^) as indicated in the proof of Proposition 10. □ 
When equality a = P holds then the simpler results follow - 

Proposition 17 For equal parameters (5 = a the probability ( [^. 76 ) and associated functions 
arising from Proposition 16 are given by 

^2(0;/) = ^2a+w ^., 2i^i(a+3/2,l;a+2;l-5^) , 

, , a+l 

S2-Fi(a + 3/2, l;a + 2; l-s2) ^ ' 

G{s) = (a + l) \s - \ 

y S2-Fi(a + 3/2, l;a + 2;l-s2j 

and for the corresponding probability ( 4- 7ll ) 

Q_ 2\2a+2 4 

a{s) = («+!){ ^2Fi(a+3/2,l;a+2;l-s2) + s 2^1(0+5/2, 1; a+2; l-s^) } ' 
G{s) = (a+2)s + 2(a+l) 



1 



s 2i^i(a+3/2, 1; a + 2; l-s2) s2Fi(a + 5/2, 1; a + 2; l-s2) 

(4.79) 

in terms of the Gauss hypergeometric function 2^1(0, b; c; z). 

Proof - we proceed in a parallel manner as in the proof of Proposition 11 and the relations 
given there, however also using the linear transformation formulae for the Hypergeometric 
functions. One such an example is 

4"S(a+l,a+l) -S2Fi(-a,i/2;3/2;s2) = i^ ^^-^'^j ^ 2i^i(Q+3/2, 1; «+2; 1-g^) . (4.80) 
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We also require the comparable contiguous relation 

-(2a + 3)s22Fi(a + 5/2,l;a + 2;l-s2) +2^1(0+3/2, l;a + 2;l-s2) +2(a + l) = . (4.81) 

□ 

One can show that these two specific cases are solutions to our third order differential 
equation ( |4.73| ), after noting the above contiguous relation linking the two Hyper geometric 
functions for the case of = 2. 

The special case of a = /? = and general N is no longer a simple case as was the 
situation for the endpoint interval. 



4.6 Jacobi to Hermite Limit for the Interior Interval 

As we have shown in Subsection 4.3 one would also expect the limit a — > 00 would recover 
the result for the Hermite ensemble on {—t, t), under the scaling in ( 4.48| ). This result is given 
in |18|, although not in a form useful for our purposes. In that reference a second order ODE 
is given for R{t) whereas we want the corresponding ODE for R{t). Using the same kind of 
elimination indicated in the proof of Proposition 4 on their Equations (5.31) and (5.32) we 
find the result 



tr" + 2r' = -2t[t - g] - 2g^ {a'f + At[t - g]a - 2Nt^[t - g]^ 



(4.82) 



where g = Vf^ + 2r' and a(t) = tr{i). To verify that ( [4.73| ) has the correct limiting behaviour, 
we first find that the substitution ( [1.48D , to leading order in o? yields 



(r'O 



tl\2 



2 2 



+ 2r'(7[t - g\ 



^ 2\ {tr"-2r')\ 



(4.83) 



and that differentiation of the second order Hermite ODE, ( 4.82| ), gives 
tr'" + 3r" = -4t + Ag'^{tr' + r- 2Nt) + 2g 

- 8t{r -Nt)g+{t + r"){2t^ + tr" + 2r')/g'^ 



(4.84) 



with g as before. By repeating the steps described in Subsection 4.3 we can show that both 
( |4.82| ) and ( [4.84 ) can be combined to give the limiting case, ( 4.83| ). 
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5 Reductions to Painleve transcendents 



The analysis in the preceding sections yielded a number of nonlinear differential equations of 
the second and third order. In this section we present their solutions in terms of the fifth and 
sixth Painleve transcendents. In all cases, we can arrive at one of the equations of the second 
order and second degree found by Chazy ||^, ^] and subsequently rederived by a number of 
authors. 

We begin with equation ( |3.19| ) for the variable R{s). This equation is equivalent un- 
der a gauge transformation to an equation first derived by Chazy by transformation of the 
Painleve-V equation. It is the fourth member of the set of five equations denoted (II) in [Q] 
and (C) in p. 342]. It was also derived by Bureau |Q, pp. 204-206] in his investigation of 
second-order second-degree equations. 

The solution of ( |3.19| ) is 

~ eiw' - 2sw N{w - 1) +ei 

R[s) = - — — , (5.1) 

Aw{w - 1) Asw 

where ei := ±1 and w{s) (not to be confused with w used previously) satisfies the differential 
equation, 

+ (!^(iV.„_(^^)+2,.._?£M!^, (5,2) 

2s"= [ W ] W — 1 

The prime denotes d/ds. Under the change of variable, x = s^, the latter equation becomes 
the standard Painleve-V equation, 

d^w f 1 11/ dw\'^ 1 dw 



+ 



dx"^ {210 w — 1 j \dx J X dx 



- 



{w-lf \ I3\ -fw 5w{w + l) . ^, 

+ ^ ^ \aw + -} + — + — ^ , 5.3 

W \ X w — 1 



with parameters, 

a = lN', p = -l(N-ei)\ 7 = 1/2^1, '5 = -1/2- (5-4) 
In terms of the same w{s), the solution of the companion equation ( |3.20D is 

8sw{w — 1)^ L J 

X ^sw' - N{w -if - {2s'^ + l)w + 1^ (5.5) 

and the auxiliary variable h{s) is given by 

2sw^-€iw' , N{w -l) + ei 
2w(w — 1) 2sw 

The last two equations also furnish the solution of ( |3.26| ) with h{s) being identified with 
the square root with the upper sign. The identity ( |3.27 ) is satisfied identically by the forms 
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of and R{s) given here. They also solve equation ( 4.49| ) after renaming the variables in 
( 4.49| ) according to t — > s, r(i) R{s). In that case, the parameter N plays the role of the 
third integration constant because ( [4 .491) can be obtained from ( 3.2C| ) by differentiating out 
the parameter N. 

On the other hand, equation ( [4.50| ), which was also obtained by differentiation of ( 3.20 ), 
is a true generalisation of ( p.20D . Its first integral is 

sR" + 2R' = 2s{s - h) 



- 2h^J{R + sR'f - 4s2(s - h)R - 2iVs2(s - hf + Ki , (5.7) 

where h denotes V s"^ — 2R' and Ki is the integration constant. Although the case of nonzero 
Ki is not relevant to the discussion in Section 4, we nevertheless have another integrable 
equation. Its solution is 



R 



^^^^^^^{W + iV(.-l)^ + (2.^-l). + l} 
X I W - N{w - if - (25^ + l)w + l} 

K{W + 1){N{W - 1)2 - 2^2^;} 

2sw{w — 1) 

^ A->-l)(3„. + l)^ (5.8) 

2sw 

, 2sw^-eiw' , {N - 2K){w - 1) + ei . ^. 

n = -. r h , 0-9) 

2w{w - 1) 2sw ' ^ ^ 

where ei := ±1, Ki = 8K'^{N — 2K), and w{s) satisfies the differential equation, 
zw w — 1 ) s 



w 



2s^w(w + 1) , , 

+ 2{AK + ei)w ' . (5.10) 

As before, the change of variable x = transforms ( 5.101 ) to a standard Painleve-V equation. 
When ii' = 0, these results reduce to the above results for ( p.20| ). By changing the sign of R, 
these formulae, with zero or nonzero Ki as appropriate, also solve equations (|4.82| )-( 4.84 ). 

We now turn our attention to equation ( 4.29| ), which presents a somewhat greater chal- 
lenge. This differential equation is of the third order and second degree. We find that it 
is more manageable when written in terms of the auxiliary variable H{s) rather than <j{s). 
After that, we observe a slight improvement by choosing a new variable y(s), in terms of 
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which we have 

H{s) = {l-s^)y-ais, (5.11) 

ais) = l{a^sW - (s' - l)y'y\sHs' - \)\yy"' - y'y") 
+ 2s{s' - 1) [(4^2 - l)yy" - {2s^ - l){y'f] 

- 2y[2s'^{s'^ - ify^ + 2,ais^{s^ - l)y - + 3s^ + l\y' 

- 2sy'^[{s^ - l)(3s2 - l)y^ + 2ais{'is^ - 2)y 

+ 2(af - l)s2]}. (5.12) 

The differential equation satisfied by y{s) takes the form, 

{2a-a,shY = -s\s^-lf{yy"-{y'f} 

-2s\s'-l)yy' + y\sHs'-lfy' 

+ 2ais^{s^ - l)y + (af - l)s^ + l}, (5.13) 

where it is understood that a is to be replaced by the right-hand side of ( ^.12| ). This 
differential equation of the third order and second degree admits the first integral, 

{.(.2-i)V + 2(.2-i)(2.2_iy 

+ 8s^y^ + Uais^y"^ + 5(25^ + Ki - 6)y|^ 

= 4{(s2 + l)y + ais}'{(s2 - l)\sy' + yf 

+ sV(4sV + 8aisy + i<Ci -4)}, (5.14) 

where Ki is the constant of integration. 

Equation ( 5.14 ) is equivalent under a gauge transformation to the fifth member of the 
aforementioned set of Chazy equations. Chazy found this equation by transforming the 
Painleve-VI equation, but he did not show the actual formula. Bureau also obtained this 
equation 1^, pp. 200-202] but did not solve it. The same equation in a different gauge was 
obtained by Fokas and Yortsos who gave the reduction to Painleve-VI. 

The variables y(s), H(s), and a{s) are given in terms of a function w{s) by the formulae, 

eis{s^ - l)w' -{w- l){ei{w + s^) + 2ais^} . ^ , 

2s(s^ — l)w 

_ eisjl - s'^)w' + €i{w - l){w + s^) - 2aig^ 

{s{s^ - l)w' - {w - l){w + s'^)}^ 
a = — — 



8sw{w — l){w — s^) 
s{w - l){{Ki - 4)w - 4a^s^} 



8w{w — s 



2\ 



(5.17) 
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where ei := ±1. The function w{s) satisfies the differential equation, 

2 [ w w — 1 w — \ 



1 2s 2s 

- + H J > w 

s — 1 w — s"^ 



^ w{'w — l){w — 



2s2(s2 _ 1)2 

(l + 2eiai)V (Ki-4«f)s2 (s2-l) 
w'^ {w — s 



Under the change of variable, x = s , the latter equation becomes the Painleve-VI equation, 



dx"^ 2 \w w — 1 w — X \ \dx J 



1 1 1 ^ dw 

- + T + 



with parameters. 



a 



X X — 1 w — X j dx 

w{w - l){w-x) [ /jrc 7(x-l) \ 
+ 27 T\2 i" + ~ + 7 T\2+7 ^f' (.^-l^) 



/?= -i(l + 2eiai)2, 7 = 0, (5 = -4Qf). (5.20) 



These formulae also solve ( [4.73| ) with the sign of a changed. 

We draw the reader's attention to another set of five second-order second-degree Chazy 
equations which make regular appearances in random matrix theory. This set is denoted (III) 
in and (B) in [|lp. 340] and has also been studied by Bureau [0, ^, Fokas and Ablowitz [^, 
Jimbo and Miwa [p^ . Appendix C], and Cosgrove and Scoufis This set can be embraced 
in a single equation, 

iy"f = 2^|ci(2;y' - yf + C2y'{xy' - yf + cz{y'f{xy' - y) 

+ CA{y'f + c^{xy' - yf + c^y' {xy - y) + C'j{y'f 

+ C8(xy'-2/) + C9y' + cio}, (5.21) 

where g{x) := cix^ + C2X^ + C3X + C4 and the prime denotes djdx. This is equation SD-I, 
introduced in pp. 57, 65-73]. It was called there the "master Painleve equation" because 
it unified all of the six Painleve transcendents into a single equation. Only four of the 
10 parameters in SD-I are essential because the equation retains its shape under gauge 
transformations of the form, 

ax + b _ hy + kx + m . . 

X = , y = ; . (5.22) 

cx + d cx + d 

By using this gauge freedom appropriately, equation SD-I can be split into six normalised 
forms, whose solutions in terms of Painleve transcendents are given in B. 
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The full version of SD-I generates an abundant supply of third-order equations under 
differentiation. For example, replacing Cj by Cj + -ftTcj+g for i = 5, . . . , 10 and differentiating 
out the parameter K produces a 15-parameter equation (after normalising cn, say, to 1 or 0). 
In a similar fashion, one can generate a 25-parameter equation of the third order and second 
degree. These big equations are all integrable in terms of Painleve-VI in the generic case and 
one of the other five Painleve transcendents or elliptic functions in the remaining cases. 

A special case of the aforementioned 15-parameter equation appears, for example, in 
Tracy and Widom [^] . This case was solved recently in terms of Painleve-VI by Haine and 
Semengue [|l^]. A more recent example, also by Tracy and Widom [^], is the equation, 

/' = l/l+^_\(y")2 _ i " 

^ 2\y'^y'-l/^^ ^ 

2{k + n) , , x + n , . 
y (y - 1) + [n-x + 2y) 

{n + yf [x-yf 



y = < w 

4(A; -|- n)w — 1 



Aw 



(5.23) 



2x'^y' 2x2(y'-l) 
the prime denoting d/dx. This equation has the first integral, 
x'^{y"f = -A{k + n)x{y'f 

+ {4{k + n)y + x'^ + 2{2k + 3n)x} {y'f 
- {2{x + 2k + 3n)y + 2nx + n^}y' 

+ {y + nf - AKiy'iy' -1), (5.24) 

where Ki is the constant of integration. 

The latter equation is gauge-equivalent to the normalised form SD-I.b in . Its solution 

is 

1 f xw' 1 ^ 

2)w \ w — 1 

{w - l){k{w - 1) + n{w + 1)} 



1 r 2nx(w - 2) x'^w , AKi ] 

4{k + n){ w-1 {w-iy w } 



where w{x) satisfies the Painleve-V equation (|5.3D with parameters, 

a = i/2(,k + n + eif, (3 = y^in^ - k^ + AK^), l = -n, <5 = -y2, (5.26) 
with ei := ±1. 
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